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ROBUSTNESS OF SCALE-FREE SPATIAL NETWORKS 


EMMANUEL JACOB AND PETER MORTERS 

Abstract: A growing family of random graphs is called robust if it retains a 
giant component after percolation with arbitrary positive retention probabil¬ 
ity. We study robustness for graphs, in which new vertices are given a spatial 
position on the d-dimensional torus and are connected to existing vertices 
with a probability favouring short spatial distances and high degrees. In this 
model of a scale-free network with clustering we can independently tune the 
power law exponent r of the degree distribution and the rate —6d at which 
the connection probability decreases with the distance of two vertices. We 
show that the network is robust if r < 2 -|- but fails to be robust if r > 3. 

In the case of one-dimensional space we also show that the network is not 
robust if T > 2 -|- This implies that robustness of a scale-free network 
depends not only on its power-law exponent but also on its clustering fea¬ 
tures. Other than the classical models of scale-free networks our model is not 
locally tree-like, and hence we need to develop novel methods for its study, 
including, for example, a surprising application of the BK-inequality. 
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1. Motivation 

Scientific, technological or social systems can often be described as complex networks of inter¬ 
acting components. Many of these networks have been empirically found to have strikingly 
similar topologies, shared features being that they are scale-free, i.e. the degree distribu¬ 
tion follows a power law, small worlds, i.e. the typical distance of nodes is logarithmic or 
doubly logarithmic in the network size, or robust, i.e. the network topology is qualitatively 
unchanged if an arbitrarily large proportion of nodes is removed from the network. Barabasi 
and Albert [2] therefore concluded fifteen years ago ‘that the development of large networks is 
governed by robust self-organizing phenomena that go beyond the particulars of the individ¬ 
ual systems.’ They suggested a model of a growing family of graphs, in which new vertices 
are added successively and connected to vertices in the existing graph with a probability 
proportional to their degree, and a few years later these features were rigorously verified in 
the work of Bollobas and Riordan, see unmiis]. 
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In the years since the publication of [2] there have been many refinements of the idea of pref¬ 
erential attachment, introducing for example tunable power law exponents jsnmaiig], node 
fitness or spatial positioning of nodes [SSKHEH]- Some of these refinements 

attempt to introduce or explain clustering, the formation of clusters of nodes with an edge 
density significantly higher than in the overall network. The phenomenon of clustering is 
present in real world networks but notably absent from most mathematical models of scale- 
free networks. The present paper investigates a spatial network model, introduced in 1211, 
defined as a growing family of graphs in which a new vertex gets a randomly allocated spatial 
position representing its individual features. This vertex then connects to every vertex in 
the existing graph independently, with a probability which is a decreasing function of the 
spatial distance of the vertices, the time, and the inverse of the degree of the vertex. The 
relevance of this spatial preferential attachment model lies in the fact that, while it is still 
a scale-free network governed by a simple rule of self-organisation, it has been shown to 
exhibit clustering. The present paper investigates the problem of robustness and is probably 
the first rigorous attempt to understand the global topological structure of a self-organised 
scale-free network model with clustering. 


In mathematical terms, we call a growing family of graphs robust if the critical parameter 
for vertex percolation is zero, which means that whenever vertices are deleted independently 
at random from the graph with a positive retention probability, a connected component 
comprising an asymptotically positive proportion of vertices remains. For several scale- 
free models, including non-spatial preferential attachment networks, it has been shown that 
robustness holds if the power law exponent satisfies r < 3, see for example mm- At a first 
glance one would maybe expect this behaviour to persist in the spatial model. It is known 
that robustness in scale-free networks relies on the presence of a hierarchically organised core 
of vertices with extremely high degrees, such that every vertex is connected to the next higher 
layer by a small number of edges, see for example m- Our analysis of the spatial model 
shows that, although a hierarchical core still exists if r < 3, whether vertices in the core are 
sufficiently close in the graph distance to the next higher layer depends critically on the speed 
at which the connection probability decreases with spatial distance, and hence depending on 
this speed robustness may hold or fail. The observation that robustness depends not only on 
the power law exponent but also on the clustering of a network appears to be new, though 
similar observations have been made for a long-range percolation model m- Unlocking this 
phenomenon is the main achievement of this paper. 


The main structural difference between the spatial and classical model of preferential attach¬ 
ment is that the former exhibits clustering. Mathematically this is measured in terms of a 
positive clustering coefficient, meaning that, starting from a randomly chosen vertex, and 
following two different edges, the probability that the two end vertices of these edges are 
connected remains positive as the graph size is growing. This implies in particular that local 
neighbourhoods of typical vertices in the spatial network do not look like trees. However, the 
main ingredient in almost every mathematical analysis of scale-free networks so far has been 
the approximation of these neighbourhoods by suitable random trees, see P [20111 [21]. As 
a result, the analysis of spatial preferential attachment models requires a range of entirely 
new methods, which allow to study the robustness of networks without relying on the local 
tree structure that turned out to be so useful in the past. Providing these new methods is 
the main technical innovation in the present work. 
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2. The spatial preferential attachment model 


While spatial preferential attachment models may be defined in a variety of metric spaces, 
we focus on homogeneous space represented by a d-dimensional torus of unit volume, given 
as Ti = (— 1 / 2 , 1 / 2 ]^^ with the metric di given by 

di(x, y) = min {d(x, y + u) \ u & 0 ) 1 }^}) for x,y € Ti, 

where d is the Euclidean distance on Let X denote a homogeneous Poisson point process 
of finite intensity A > 0 on Ti x (0, oo). A point x = (x, s) in A is a vertex x, born at time s 
and placed at position x. Observe that, almost surely, two points of X neither have the same 
birth time nor the same position. We say that (x, s) is older than (y, t) if s < t. For t > 0, 
write Xf for X n (Ti x (0,t]), the set of vertices already born at time t. 

We construct a growing sequence of graphs {Gt)t>o, starting from the empty graph, and 
adding successively the vertices in X when they are born, so that the vertex set of Gt 
equals Xt. Given the graph Gt- at the time of birth of a vertex y = (y, t), we connect y, 
independently of everything else, to each vertex x = (x, s) G Gt-, with probability 

where Z{x,t—) is the indegree of vertex x, defined as the total number of edges between x 
and younger vertices, at time t—. The model parameters in Q are the attachment rule 
/: N U {0} —)• (0,oo), which is a nondecreasing function regulating the strength of the 
preferential attachment, and the profile function (p: [ 0 , oo) —>■ ( 0 , 1 ), which is an integrable 
nonincreasing function regulating the decay of the connection probability in terms of the 
interpoint distance. 

The connection probabilities in (j^ may look arcane at a first glance, but are in fact com¬ 
pletely natural. To ensure that the probability of a new vertex connecting to its nearest 
neighbour does not degenerate, as t t oo, it is necessary to scale di(x,y) by which is 

the order of the distance of a point to its nearest neighbour at time t. The linear dependence 
of the argument of <p on time ensures that the expected number of edges connecting a new 
vertex to vertices of bounded degree remains bounded from zero and infinity, as t f oo, as 
long as if is integrable on [ 0 ,oo), or equivalently x i—(^(d(x, 0 )'^) is integrable on M'^. 

The model parameters A, / and g? are not independent. Indeed, if f (p(d(x, 0)'^) dx = y, > 0, 
we can modify (f to (p o (yid) and / to yf, so that the connection probabilities remain 
unchanged and 

J ip{d{x,0)'^) dx = 1 . (2) 

Similarly, if the intensity of the Poisson point process T is A > 0, we can replace X by 
{(x. As): (x, s) G T} and / by A/, so that again the connection probabilities are unchanged 
and we get a Poisson point process of unit intensity. From now on we will assume that both 
of these normalisation conventions are in place. 


Under these assumptions the regime for the attachment rule / which leads to power law 
degree distributions is characterised by asymptotic linearity, i.e. 


lim 


m 

k 


= 7, 


for some 7 > 0. We henceforth assume asymptotic linearity with the additional constraint 
that 7 < 1 , which excludes degenerate cases with infinite mean degrees. 
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We finally assume that the profile function (p is either regularly varying at infinity with index 
— 5, for some 5 > 1, or decays quicker than any regularly varying function. In the latter 
case we set 6 = oo. Intuitively, the bigger 5, the stronger the clustering in the network. Our 
assumptions, in particular the assumption that p does not take the values 0 or 1 , help us 
avoid some geometric constraints that are not of major interest. See Figure 1 for simulations 
of the spatial preferential attachment network indicative of the parameter dependence. 



Figure 1. Four simulations of the network in the case d = 2, based on the 
same realisation of the Poisson process, with parameters (clockwise from top 
left) (a) 7 = 0.5 and 6 = 2.5, (b) 7 = 0.75 and 6 = 2.5, (c) 7 = 0.5 and <5 = 5, 

(d) 7 = 0.75 and <5 = 5. The pictures zoom into a typical part of the torus. 

A similar spatial preferential attachment model was introduced in [T] and studied further 
in [28l [13]. In this model it is assumed that the profile functions has bounded support, 
more precisely p = for P ^ (0,1] and r satisfying ([^. This choice of prohle function, 

roughly corresponding to the boundary case d t c>o, is too restrictive for the problems we 
study in this paper, as it turns out that robustness does not hold for any value of r. There 
are also spatial long-range percolation models which have a qualitatively similar behaviour 
to our networks, see for example [mini [26], but these models are easier to analyse and the 
methods of this analysis are quite different. 

Local properties of the spatial preferential attachment model were studied in [27] , where this 
model was first introduced. It is shown there, among other things, that 

• The empirical degree distribution of Gt converges in probability to a deterministic 
limit p. The probability measure ^ on {0} U N satisfies 

pi{k) = as k 00 . 

In other words, the network {Gt)t>o is scale-free with power-law exponent r = 1 
which can be tuned to take any value r > 2 . See m Theorem I and 2]. 
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• The average over all vertices v € Gt of the empirical local clustering coefficient at v, 
defined as the proportion of pairs of neighbonrs of v which are themselves connected 
by an edge in Gt, converges in probability to a positive constant > 0, called the 
average clustering coefficient. In other words the network {Gt)t>o exhibits clustering. 
See [271 Theorem 3]. 


3. Statement of the main results 

We now address the problem of robustness of the network {Gt)t>o under percolation. Recall 
that the number of vertices of the graphs Gt, t > 0, form a Poisson process of unit intensity, 
and is therefore almost surely equivalent to t as t\ oo. Let Gt C Gt be the largest connected 
component in Gt and denote by \Gt\ its size. We say that the network has a giant component 
if Gt is of linear size or, more precisely, if 


lim lim sup P 
sfO t^oo 



= 0 . 


We say it has no giant eomponent if Gt has sublinear size or, more precisely, if 


lim inf P ( < e 

t^ao y t 


1 for any e > 0. 


If G is a graph with vertex set X, and p G (0,1), we write for the random subgraph of G 
obtained by Bernoulli percolation with retention parameter p on the vertices of G. We also 
use 'PX for set of vertices surviving percolation. The network {Gt)t>o is said to be robust if, 
for any fixed p G (0,1], the network {PGt)t>o has a giant component and non-robust if there 
exists p G (0,1] so that {PGt)t>o has no giant component. 


Our main result concern phases of robustness or non-robustness for the spatial preferential 
attachment network. In classical non-spatial preferential attachment models, there is a phase 
transition for robustness when r crosses the critical value 3, see [21]. It is easy to believe that 
the spatial structure does not help robustness. Our main result shows that in the spatial 
model robustness is still possible, but at least in the case d = 1 the phase transition does 
not occur at r = 3, but at a smaller value depending on 6. 


Theorem 1. The spatial preferential attachment network {Gt)t>o is 

(a) robust if ■y > or, equivalently, if t <2 + ^; 

(b) non-robust if j < ^ or, equivalently, if t > ?>. 

In the case of one space dimension, d = 1, the network is 

(c) non-robust if 'j < or, equivalently, if t > 2 + . 

Remark 2. For a suitable range of parameters, this seems to he the first instance of a seale- 
free network model which combines robustness with clustering features. We conjecture that 
nonrobustness occurs in any dimension if j and thus the critical value for r equals 

2 -|- y, but our proof techniques do not allow to prove this, see Figure 2 for a phase diagram. 
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Figure 2. The densely shaded area corresponds to the known robustness 
phase 7 > (5/((5 + 1). The lightly shaded area corresponds to the known non¬ 
robustness phase. In between, no result of either robustness or non-robustness 
is known, though non-robustness is conjectured. 


Remark 3. Our approach also provides heuristics indicating that in the robust phase 5 {t — 
2) < 1 the typical distances in the robust giant component are asymtotically 


(4 + 0(1)) 


log log t 

-log(5(r-2))’ 


namely doubly logarithmic, just as in some nonspatial preferrential attachment models. The 
constant coincides with that of the nonspatial models in the limiting case d 1, see [23l ITS] . 
and goes to infinity as 6{t — 2) —)• 1. It is an interesting open problem to confirm these 
heuristics rigorously. 


4. The limit model and proof strategies 


Before describing the strategies of our proofs, we briefly summarise the techniques developed 
in m in order to describe the local neighbourhoods of typical vertices by a limit model. We 
will heavily rely on these techniques in the present paper. 


Canonical representation. We first describe a canonical representation of our network {Gt)t>o- 
To this end, let T be a Poisson process of unit intensity on Ti x (0, oo), and endow the point 
process X x X with independent marks which are uniformly distributed on [0,1]. We denote 
these marks by Vx,y or V(x,y), for x,y G T. 


If T C Ti X (0, oo) is a finite set and W: T x T —)• [0,1] a map, we define a graph G^{y, W) 
with vertex set y by establishing edges in order of age of the younger endvertex. An edge 
between x = {x,t) and y = {y, s), t < s, is present if and only if 


where Z(x,s—} is the indegree of x at time s—. A realization of X and V then gives rise to 
the family of graphs {Gt)t>o with vertex sets Xt = Xn (Ti x ( 0 , t]), given by Gt = G^{Xt, V), 
which has the distribution of the spatial preferential attachment network. 


Space-time rescaling. The construction above can be generalised in a straightforward manner 
from Ti to the torus of volume t, namely T* = (—^ t^^^, \ equipped with its canonical 

torus metric dj. The resulting functional, mapping a finite subset T C T* x (0,oo) and a 
map from T x T —^ [ 0 , 1 ] onto a graph, is now denoted by G*. 
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We introduce the rescaling mapping 

hf. Tix(0,t] ^ TiX(0,l], 

{x,s) l-A 

which expands the space by a factor the time by a factor 1/t. The mapping ht operates 
on the set X, but also on V, by ht(V)/ij(x),ht(y) := ^x,y The operation of ht preserves the 
rule Q, and it is therefore simple to verify that we have 

G\ht{Xt),ht{V)) = ht{G\Xt,V)) = ht{Gt), 

that is, it is the same to construct the graph and then rescale the picture, or to first rescale 
the picture, then construct the graph on this rescaled picture. Observe also that ht{Xt) 
is a Poisson point process of intensity 1 on x (0,1], while htiV) are independent marks 
attached to the points of ht{Xt) x ht{Xt) which are uniformly distributed on [0,1]. 


Convergence to the limit model. We now denote by T a Poisson point process with unit 
intensity on x (0,1], and endow the points oi X xX with independent marks V, which are 
uniformly distributed on [0,1]. For each t > 0, identify (—^ ^ t^/^/2Y and Tj, and write 

X^ for the restriction of Al to x (0,1], and V* for the restriction of V to T* x AfF In the 
following, we write G^ or G*(A’,V) for G*(A’*,V*). We have seen that for fixed t G (0,oo), 
the graphs G* and ht{Gt) have the same law. Thus any results of robustness we prove for 
the network (G*)t>o also hold for the network (Gt)i>o. It was shown in [271 Proposition 5] 
that, almost surely, the graphs G* converge to a locally finite graph G°° = G“(A’,V), in 
the sense that the neighbours of any given vertex x G Af coincide in G* and in G“, if t is 
large enough. It is important to note the fundamentally different behaviour of the processes 
iG%>o and {Gt)t>o- While in the former the degree of any fixed vertex stabilizes, in the 
latter the degree of any fixed vertex goes to oo, as 11 oo. We will exploit the convergence of 
G* to G°° in order to decide the robustness of the finite graphs G*, and ultimately Gt, from 
properties of the limit model G°°. 


Law of large numbers. We now state a limit theorem for the graphs ^G* centred in a randomly 
chosen point. To this end we denote by the law of AfjV together with independent 
Bernoulli percolation with retention parameter p on the points of X. For any x G x (0,1] 
we denote by ^Px the Palm measure, i.e. the law ^P conditioned on the event {x G ^X}. 
Note that by elementary properties of the Poisson process this conditioning simply adds the 
point X to ^X and independent marks Vx,y and Vy^x, for all y G A", to V. We also write ^Ex 
for the expectation under PPx. Let ^ ^ (x, G) be a bounded functional of a locally-finite 

graph G with vertices in x (0,1] and a vertex x G G, which is invariant under translations 
of W^. Also, let (x, G) be a bounded family of functionals of a graph G with vertices 

in Tt X (0,1] and a vertex x G G, invariant under translations of the torus. We assume 
that, for U an independent uniform random variable on (0,1], we have that ^*((0, G), ^G*) 
converges to f,{{0,U),PG°°) in PP(o,! 7 )"Probability. Then, in ^’P-probability, 


] E «< "<?<) 

Vt 


■ D 


PTE , 


(4) 


This law of large numbers is a minor modification of the one given in m Theorem 7], which 
covers the case p = 1. 
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4.1. Robustness: strategy of proof. 

Existence of an infinite component in the limit model. We first show that, under the assump¬ 
tions that 7 > and p G ( 0 , 1 ], the percolated limit model has an infinite connected 
component. This uses the established strategy of the hierarchical core. The young vertices, 
born after time are called connectors. Fix a G (1, Starting from a sufficiently old 

vertex xq G pG°°, we establish an infinite chain {:x.k)k>i of vertices = {xk,Sk) such that 
Sk < 'Sfc-i) be. we move to increasingly older vertices, and Xfc_i and x^ are connected by a 
path of length two, using a connector as a stepping stone. 

Transfer to finite graphs using the law of large numbers. To infer robustness of the network 
iG%>o from the behaviour of the limit model we use Q on the functional ^i(x, G) defined 
as the indicator of the event that there is a path in G connecting x to the oldest vertex 
of G. We denote by ^(x, G) the indicator of the event that the connected component of x is 
infinite and let 

^9 := [ (0 ,^)1 the component of (0, u) in PG°° is inhnite } du. (5) 

JO 

If lim^t((0, U),PG^) = ^((0, U), PG°°) in probability, then the law of large numbers 0 implies 

^ ^fi^,PG^)^pP9. 

The sum on the left is the number of vertices in pG^ connected to the oldest vertex, and we 
infer that this number grows linearly in t so that a giant component exists in (^G*)t>o. This 
implies that (G*)t>o and hence {Gt)t>o is a robust network. However, while it is easy to see 
that limsup^.|.oo ^*((0; U),PG^) < ^((0, U), pG°°), checking that 

lirninf 6 (( 0 ,C/),^'G*) > ^(( 0 , H), ( 6 ) 

tfoo 

is the difficult part of the argument. 

The geometric argument. The proof of Q is the most technical part of the paper. We first 
look at the hnite graph ^G* and establish the existence of a core of old and well-connected 
vertices, which includes the oldest vertex. Any pair of vertices in the core are connected by 
a path with a bounded number of edges, in particular all vertices of the core are in the same 
connected component. This part of the argument is similar to the construction in the limit 
model. We then use a simple continuity argument to establish that if the vertex (0, U) is in 
an infinite component in the limit model, then it is also in an inhnite component for the limit 
model based on a Poisson process X with a slightly reduced intensity. In the main step we 
show that under this assumption the vertex ( 0 , C/) is connected in ^G* with reduced intensity 
to a moderately old vertex. In this step we have to rule out explicitly the possibilities that 
the inhnite component of PG°° either avoids the set of eligible moderately old vertices, or 
connects to them only by a path which moves very far away from the origin. The latter 
argument requires good control over the length of edges in the component of (0, U) in pG°°. 
Once the main step is established, we can hnally use the still unused vertices, which form a 
Poisson process with small but positive intensity, to connect the moderately old vertex we 
have found to the core by means of a classical sprinkling argument. 
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In Section 5 we will carry out this programme and prove robustness. In fact, we sha 
that under our hypothesis 7 > a stronger statement holds, see Proposition 


14 


1 prove 
which 

also implies that the size of the second largest connected component in ^Gt does not grow 
linearly. In other words, we will see that in this regime the network has a unique giant 
component. 


4.2. Non-robustness: strategy of proof. 

Using the limit model. If 7 < ^ it is very plausible that the spatial preferential attachment 
network is non-robust, as the classical models with the same power-law exponents are non- 
robust PH] and it is difficult to see how the spatial structure could help robustness. We 
have not been able to use this argument for a proof, though, as our model cannot be easily 
dominated by a non-spatial model with the same power-law exponent. Instead we use a 
direct approach, which turns out to yield non-robustness also in some cases where 7 > j- 
The key is again the use of the limit model, and in particular the law of large numbers. We 
apply this now to the functionals defined as the indicator of the event that the 

connected component of x has no more than k vertices. Clearly, lim^'i-oo = 

^^'“'((0, U), PG°°), and therefore 

The left hand side is asymptotically equal to the proportion of vertices in ^G^ which are in 
components no bigger than k. As fc —)> 00 the right hand side converges to 1 — ^9. Hence 
if = 0 for some p > 0, then (G*)t>o and hence {Gt)t>o is non-robust. It is therefore 
sufficient to show that, for some sufficiently small p > 0, there is no infinite component in 
the percolated limit model 


Positive correlation between edges. We first explain why a naive first moment calculation 
fails. If (0, U) has positive probability of belonging to an infinite component of pG°° 
then, with positive probability, we could find an infinite self-avoiding path in pG°° start¬ 
ing from xq = (0,t/). A direct first moment calculation would require to give a bound 
on the probability of the event {xq-h-xi-h- • • • -H-Xn} that a sequence (xq, ..., x^) of distinct 
points Xj = (xj. Si) conditioned to be in X forms a path in G°°. If this estimate allows us to 
bound the expected number of paths of length n in G°° starting in xq = (0, U) by C'"', for 
some constant C, we can infer with Borel-Cantelli that, if p < l/C, almost surely there is 
no arbitrarily long self-avoiding paths in PG°°. For a variety of non-spatial models the event 
{xQ-H-xi-f-)- • • • -fAx^} can be decomposed into independent, or negatively correlated, events of 
the form {xj-fAXj+i}, or {xj-i-fAXiO-Xj+i} with Si < Si-i, Si, the probability of which can be 
easily estimated, see for example [T8| . For spatial networks however such a decomposition is 
not possible. Indeed, the events {xjOXj+i} and {xfc-H-Xfc+i} are not independent if the in¬ 
terval / = (sj, Sj+i)n(sfc, Sfc+i) is nonempty, because the existence of a vertex in A’n(M‘^x/) 
which is relatively close to both xj and Xk is likely to connect to both of these vertices and 
make their indegrees grow simultaneously. Observing that all the events {x^ox^+i} are 
increasing in the Poisson point process X, we can argue by Harris’ inequality that they are 
positively correlated. Because the positive correlations play against us, it seems impossible 
to give an effective upper bound on the probability of a long sequence to be a path, therefore 
making this first moment calculation impossible. 
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Quick paths, disjoint occurrence, and the BK inequality. As a solution to this problem we 
develop the concept of quick paths. Starting from a sequence (xo,...,x„) in with 

xq = X and x„ = y we construct a new sequence (zq, ..., z^), with zq = x and z^, = y, such 
that at least half of the points are in pG°°, and the remaining ones are in G°°. This sequence, 
called a quick path, has the property that the event {zQ-f-)- • • • -H-Zm} can be split into smaller 
parts, in the sense that it implies the disjoint occurrence of events {zj-fA • • • •fAZj+ 4 } involving 
five or fewer consecutive vertices of the sequence. The concept of disjoint occurrence is due 
to van den Berg and Kesten, and the famous BK-inequality states that the probability of 
events occurring disjointly is bounded by the probability of their product. It is tedious, but 
not hard, to estimate the probability of these paths of length no more than five, and the 
estimate produces the necessary bounds to complete the argument. 

Instead of defining quick paths and disjoint occurrence here, we just give a flavour by showing 
how to deal with a path xq 0 x 10 x 2 , for xo,xi,X 2 G pG°° and so,S 2 < si. To move to the 
quick path we let zq = xq, Z 2 = X 2 and replace the vertex xi G pg°° by the oldest vertex 
zi = {zi,ui) G G°° such that sq, S 2 < ui and Z 0 OZ 1 OZ 2 . Now any vertex z[ = G G°° 

with So, S 2 < u'l < can only influence the indegree of either zq or Z 2 at time ui, but never 
both. This means, loosely speaking, that zoOziOZ 2 being a quick path implies the disjoint 
occurrence of the events {zqOzi} and {Z 1 OZ 2 }. 

In Section 6 we will carry out this programme and prove non-robustness if 7 < |, or if d = 1 
and 7 < . Some auxiliary lemmas used in various parts of our proofs have been postponed 

to an appendix, see Section 7. 

Summary of standard notation. We use the Landau symbols o{t), 0{t), and n(t). If (A(t))i>o 
is a family of events, we say A{t) holds with high probability, or whp{t), if the probability of 
A{t) goes to I, as 11 00 . We say A{t) holds with extreme probability, or wep{t), if it holds 
with probability at least 1 —exp(—n(log^ t)), as 11 oo. When the parameter is clear, we write 
whp or wep for whp{t) or wep{t). Observe that, if {A{t)n)n>o is a sequence of events that 
simultaneously hold wep{t), in the sense that sup„ Prob(A(t)jj) = exp(—0(log^ t)), as 11 00 , 
and n{t) satisfies n{f) = o{t^) for some x > 0, then nfc<n(t) bolds wep{t). Informally, 

the intersection of a polynomial number of events, which each holds wep, holds wep. 

5. Proof of robustness 

In the following three subsections we study percolation on the infinite graph G°°. The giant 
component for the sequence {G^)t>o of finite graphs is studied in Subsection 5.4. 

5.1. Infinite component of the infinite graph. In this section we prove that the infinite 
graphs PG°° cannot contain more than one infinite component. 

Proposition 4. In the graphs ^G°°, for p G (0,1], the number of infinite components is 
always either almost surely equal to zero, or almost surely equal to one. 

The analogue of this proposition for percolation on the integer lattice, or the Poisson random 
connection model in is known for some time. Our proof follows the classical technique 
of Burton and Keane, see |12) . We focus on the case p = 1, the other cases being similar. 
The first step of this proof is to use the ergodicity of the model to deduce that the number 
of infinite components is an almost sure constant in N U {0, 00 }. The second step is to 
ensure that this constant cannot be a finite number k >2. Informally, if the graph contains 
k > 2 infinite components and if x and y are two vertices belonging to 2 different infinite 
components, then, sampling again the random variable Vx,y gives us a positive probability 
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of connecting these two components, hence decreasing the number of infinite components by 
(at least) one, leading to a contradiction. 

We focus on the last step, which is to ensure that the number of connected components 
cannot be infinity. We suppose, for contradiction, that 

{HI) almost surely, the graph G°° contains infinitely many infinite components. 


We say a vertex x G is a trifurcation if it is linked to at least three other vertices 
xi,X 2 ,X 3 G G°°, so that if x (and all its adjacent edges) is removed, then xi, X 2 , X 3 are in 
three different infinite components of G°°. Note that prior to the removal of x the vertices 
xi, X 2 , X 3 are all in the same infinite component as they are all connected to x. 

Lemma 5. If (HI) is satisfied, then / dnP(o,u) ((0) ^) ® trifurcation^ > 0. 

Jo 


Proof. We write P* {du dco) = du P{o,ii) {doj) for the underlying measure and note that 
U = U{u,u}) = u is the uniformly distributed birth time of the vertex located at the origin. 
Recall that G°° = G°°{X,V) and abbreviate := G°°{X — {(0,r)}, V), so that the law of 
under P* is the same as the law of G°° under P. Observe that the conditional probability 
given Gg° and U, that the vertex located at the origin has degree 0 in G°°, is almost surely 
in (0,1). This observation uses in particular the fact that (p does not take the values zero 
or one. Similarly, the conditional probability that the neighbouring vertices of (0, U) in 
are exactly some given xi, X2, X3 in Gg°, is also almost surely in (0,1). 

Assuming {HI), the graph contains almost surely infinitely many infinite components, 
and we may then specify arbitrarily three vertices xi, X 2 , X 3 belonging to three different ones. 
Then, under P*, there is positive probability that (0, [/) is older than these three vertices, 
and connected by an edge to exactly these three vertices. If this happens, then the presence 
of the new vertex ( 0 , rt) and the three edges linking to xi, X 2 and X 3 does not change the 
indegree of any of the other vertices. In other words, it does not interfere with the rest of 
the graph. The graphs and G°° are exactly the same except that the latter contains one 
vertex and three edges more. Therefore (0, tt) is a trifurcation of □ 


By stationarity the expected number of trifurcations in the ball B{Q, t), centered at the origin 
and of radius t, is proportional to the volume of the ball, that is to Let Et be the set of 
edges connecting a vertex x inside to a vertex y outside of B{0,t). 


Lemma 6. The cardinality of Et exceeds the number of trifurcations inB{0,t) by at least 2. 


Proof. To each trifurcation x = (x, s) with x G B{0,t), we can associate a partition of Et 
into three sets such that two edges in Et that are not in the same set are not in the same 
component of the graph with x and its incident edges removed. We obtain a compatible 
collection of partitions in the sense of Burton-Keane and the result follows accordingly. □ 

By the lemma, the expected number of elements in Et is at least proportional to But 
the expected degree of a given vertex (with birth time uniform in ( 0 , 1 ]) is finite, and so the 
expected number of neighbours at distance at least L is decreasing to zero, as L —)• 00 . It 
follows that the expectation of Et must be o{t‘^), and we obtain a contradiction. 


5.2. Continuity of the density of the infinite component. In this subsection, we are 
interested in the continuity properties of ^9 , as defined in ([^ , with respect to the parameters 
of the model. To this end we now suppose that provides a consistent family of Poisson 
point processes (aA’)o<a<i of intensity A so that, for X < p, we have \X C and hence 
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G*{\X, V) is a subgraph of G^{^X, V). For any x G x (0,1] we denote by the law 
conditioned on the event {x Denoting \G* := G^{xX, V) we again obtain convergence 

to a limit graph xG°°. We let ^G°° be the percolated limit graph. Recalling that ^(x, G) is 
the indicator of the event that the connected component of x in G is infinite, we define 

:= [' lE(^o,u)mO,n),p^)]du. 

Jo 

Proposition 7. For fixed p < 1, the function A i— )• ^0 is non-decreasing, right-continuous, 
and left-continuous everywhere except possibly at Xfip) ■= sup{A: ^6 = 0}. 

Remark 8. A similar result holds for the function p i—?• ^0, and is a variant of well-known 
results of percolation theory. 

Proof of Proposition^ We only prove left-continuity at A > Ac in the case p = 1, as the 
other parts of the statement will not be used in the sequel. Fix Ac < p < A, and let G 
be the infinite component of ^G°° . On the event that (0, u) is in the infinite component of 
xG°^, it is connected in this graph to a vertex of G. Thus there exists k € N such that the 
fc-neighbourhood of (0, u) in the graph xG°°, i.e. the set of vertices of the graph within graph 
distance k of (0, n), intersects G. Given k, with probability one, the ^-neighbourhood of 
(0, u) is the same in xG°° and in x'G°° for A' close enough to A. Hence ^((0, u), G°°{x'X, V)) 
converges almost surely, and hence in probability, to ^((0, n), G°^{xX, V)) when A' f This 
proves the left-continuity at A > Ac. □ 


5.3. Robust percolation in the limit model. In this subsection, we work in the super¬ 
critical phase 7 > We want to show that the infinite graph G°°, as well as its thinned 
versions pG°°, for any p > 0, do contain an infinite component. 

We sketch a simple strategy to get an infinite path in G°°. Start from a sufficiently old 
vertex. Use a vertex born after time | as a stepping stone to connect the old vertex by 
two edges to a much older vertex. Keep going forever, moving to older and older vertices. 
To ensure that this procedure generates an infinite path with positive probability we need 
to show that an old vertex is wep at graph distance two in G°° from a much older vertex, 
so that the failure probabilities sum to a value strictly less than one. To get the necessary 
estimates we also need to avoid using old vertices that have an exceptionally small degree. 
The expected degree of a vertex with birth time s is of order and the lemma below 
allows to define a notion of good vertices in such a way that (i) every good vertex satisfies a 
lower bound on the degree at time and (ii) most old vertices are good. To this end fix a 
function g as in Lemma of the appendix, we will never need to know more about it than 
the fact it grows slower than polynomially. 


Definition 9. For a vertex x = (x, s) in G°° we denote by Z(x, t) its indegree at time t > s, 
or in other words the number of vertices born during the time interval [s, t] that are connected 
by an edge to x. The vertex x = (x, s) is a good vertex if s < 1/2 and 

Z(x, 1/2) > s~'^/g{s~^). 


Lemma 


24 


ensures that under the vertex x with birth time s < 


2 is good whp{l/s). 


Definition 10. A vertex x = (x,s) born at time s <1/2 is locally good if its indegree in 
the graph G°°{X n ([x — x (0,1/2]) , V) is at least equal to s~'^/g{s~^). 
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The advantage of this more restrictive dehnition is that we can ensure that a vertex is locally 
good by only watching for the set of vertices nearby, up to distance Moreover, by 

Lemma 24, a vertex (x,s) is still locally good whp{l/s). The next lemma quantifies how 


young vertices allow to connect old vertices that have reached a high degree. 


Lemma 11 (Two-connection lemma). Let xi and X2 be two vertices of born before 
time 1/2. Write I := d(xi,X2)'^ and zi := Z(xj,l/2) assuming I = Ll{zi), for i G {1,2}. 
Suppose there exists e > 0 such that 

ziz^W^ = . 

Conditional on the restriction ofG°° to vertices born before time 1/2, wep(zi), the vertices xi 
and X2 are eonnected through a vertex born after time 1/2, in the graph G°°. The analogous 
result holds also for the thinned graphs pG°°, for any given p > 0. 


Proof. From the construction rule and the hypothesis Z(xi, 1/2) = zi, it is clear that every 
vertex x = (x,s) of T satisfying the conditions s > 1/2, d(x, xi)'^ < f(zi) and Vx,xi < <y?(l) 
is connected to xi. The number of such vertices is a Poisson variable with parameter of 
order zi, therefore it is wep{zi) of order zi. The probability that x connects to X2 is 


/ sd(x,X2)-^ \ 
\f{Z{x2,S-))) 


> T 


f(z2) 


On the right hand side, the numerator of the argument of is 0{i), and the Potter 
bounds for regularly varying functions, see [HI Theorem 1.5.6], ensure the right hand side 
is for any e' > 0. Hence the event that one of the vertices x is connected 

to X2 by an edge is stochastically bounded from below by a binomial random variable with 
parameters of order ^(zi) (number of trials) and Tl{{£zf^)~^~^ ) (success probability). If e' 
is chosen small enough, the expectation of this random variable is Ll{zi‘^), and therefore it 
is wep{zi) of order Ll{z]^'^). In particular, it is positive, and it stays wep{zi) positive after 
percolation with any retention parameter p > 0. □ 


Corollary 12. Suppose 7 > Choose 

first a e {l,j^^), then 13 ^ {a,1{l + a6)). 

If the vertices x = (x, s) and y = {y, t) are good vertices with t < s^^ and d(x, y)'^ < 
then they are conneeted through a vertex born after time 1/2, wep(l/s). 


With this corollary in hand, we can state and prove the following key proposition. 

Proposition 13 (Chain of ancestors). If xq = (xo,so) is a locally good vertex, then 

(1) wep(l/so) there exists a loeally good vertex xi = (xi,si) G T with si < Sq and 
d(xo,xi)'^ < Sq^. We say xi is an ancestor o/xq. 

(2) wep(l/so) there exists an infinite chain of ancestors (xfc)fc>i, namely loeally good 

vertiees satisfying and d{xk,Xk+iy < s//^ for every k>0. 

(3) wep(l/so), two consecutive aneestors of the infinite chain of ancestors are always 
eonnected through a vertex born after time 1/2, and therefore within graph distance 
two in G°°. 
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Proof. The only difficult part is to explain how to hnd the ancestors. We have ensured that 
xq is locally good, by looking only at the vertices in the set 

n {[xo - xo + X (0,1/2]). 

We always search for ancestors by moving to the right on the first coordinate. Let £ = (3 —a. 
Take /Q\ — 1 disjoint intervals of length 6 sq“^‘* inside Write for the 

centre of the fc-th interval and := (a^, 0,..., 0) + [—The blocks x + 
are disjoint and have not been observed so far. Therefore, independently of everything else, 
the probability that X (P{x + Ak) contains a locally good vertex at distance less than 
from the centre of the block with birth time in (sq/ 2 , Sg) is bounded from zero, say by c > 0 . 
One of the L'So^^'^/ 6 j — 1 independent trials with success probability > c has to succeed, 
w;ep(l/so), which proves (1). Similarly, given the k first ancestors, we hnd the {k + l)-th 
ancestor wep{l/sk). Therefore we easily see that we can hnd an inhnite chain of ancestors, 
wep{l/so), which proves (2) as well. □ 

It follows directly from Proposition that the inhnite graph percolates in the supercritical 
phase. The same proof also holds for the thinned inhnite graphs so that they also 

percolate. This immediately implies that ^6 > 0, for any p > 0. Of course, the only inhnite 
component of PG°° is a subgraph of that of G°°. In this sense, the inhnite component of the 
inhnite graph exists, is unique and robust. 


5.4. Robustness of the giant component. We now show the following result. 

Proposition 14. Let 7 > and p > 0. With high probability, the largest component of 
PG^ contains {pQ + o(l))f vertices, while the second largest contains only o{t) vertices. Henee 
there is a unique giant component, which has asymptotic density pQ > D. 

It is tempting to believe that the result follows from the robustness of the inhnite component 
of the inhnite graph by a pure approximation argument. However, the equivalence of the 
existence of a unique inhnite cluster in the inhnite graph, and that of a unique giant compo¬ 
nent in the hnite graphs, does not always hold for long-range percolation models. In the rest 
of this subsection we prove this equivalence with a signihcant effort, and our argument uses 
specihc features of the graphs in the phase 7 > The proof is organized in three parts. 
First, a direct study of the hnite picture shows that the old vertices form a core of very-well 
connected vertices. The method is similar to that of the ‘chain of ancestors’ argument for 
the limit model. Then, we explain how the law of large numbers allows to pass from the 
inhnite to the hnite picture. Finally, we have to show a convergence result for a carefully 
chosen functional to complete the proof. 


Core of well-conneeted vertices. With high probability, the oldest vertex of G^ has birth time 
within (^ 7 ^, Then Lemma ensures it is also good with high probabiht}0 We now 

work implicitly on this event. In particular, we say that a statement holds wep while the 
precise formulation should be ‘wep, the statement holds or the oldest vertex is not a good 
vertex with birth time within (^ 7 ^, ^^)’. Such an abuse of notation is used in the following 
proposition. For its formulation define, for k > 0, the 2k-eore to be the set of good vertices x 
with birth time s < t~^^°‘ . 


^The reader may note that the function log used here and in Lemma 24 plays no special role other than 


being a function growing to infinity slower than polynomially. And indeed, here and in Lemma [M} we could 
have replaced it by any other function growing to infinity slower than polynomially. 
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Proposition 15. 

(1) In G^, wep, every vertex of the 2k-core is connected, through a vertex born after 
time 1/2, to an ancestor or to the oldest vertex (or is the oldest vertex). 

(2) In G*, wep, every vertex that has reached degree at least at time 1/2 is connected 
through a vertex born after time 1/2 to the 2k-core. 


Proof. For item (1), proceed as in Proposition 13 to ensure that, wep{l/s), a good vertex 
(x, s) with birth time s < has an ancestor and is connected to it through a vertex 

born after time 1/2. The slight difference is that we work here in the finite graphs and on 
the torus. The proof has to be adapted when s~^^^ > ^ because then the blocks we 

consider cover the whole torus and overlap. In that case, recall that the oldest vertex is a 
good vertex, within distance ^ of x, and born before time Corollary 1 1 2| then ensures 
it is wep connected, through a vertex born after time 1/2, to (x, s). Finally, item (2) follows 
easily from a further application of Lemma 11 □ 


Proposition 15 has an important consequence, namely, wep, all the vertices of the 2A:-core, 
as well as all the vertices that have reached at least degree at time 1/2, belong to 

the same connected component of G^. Moreover, any two such vertices are within distance 
4:k + 4. The connected component of the core is a natural candidate for the giant component 
of the graph G*. 


Use of the laws of large numbers. Recall that Formula Q gives, whp, the asymptotic density 
of vertices of ^G^ belonging to components of size < k, for any given k £ N. As A; t oo it goes 
to 1 — ^0 and hence, whp an asymptotic proportion {1 — ^9) of vertices belongs to finite-size 
components. The remaining PQt vertices belong to large clusters, whose sizes grow with t. 
However, at this point, nothing guarantees that they form one giant component. They could 
belong to various components of logarithmic size, for example. 


To see why this is not the case, we search for an indicator function taking the value one 
on exactly one component of ^G^, which converges in probability to the indicator function 
of the event that x belongs to the infinite component of ^G°°. Inspired by the description 
of the core, we define to be the indicator of the event that x is connected through a path 
to the oldest vertex of If we can prove the convergence in probability of ^*((0, U), ^G*) 
to ^((0, U), PG°°), then the law of large numbers gives ^ ^G*) -£■ ^9. In other 

words, whp the component of the oldest vertex contains {^9 + o(l))t vertices, as claimed. 
By the previous paragraph the second largest component cannot contain an asymptotically 
positive proportion of all the vertices. This completes the proof of Proposition |14| subject to 
the assumed convergence, which we now prove. In the proof we assume p = 1 to lighten the 
notation; the general case follows by the same line of arguments. 

Convergence of to f. Recall the notation W‘^:{duduj) = duP(o^u)(da;) and U{u,uj) = u. 
We have to prove the following two statements: 

P*(At\Aoo) —^ 0 and P*(Aoo\At) —^ 0, 

t^OO t^OO 

where At denotes the event that (0, U) is connected to the oldest vertex of G*, and A^o denotes 
the event that (0,17) belongs to the infinite component of G°°. The first statement follows 
directly from the almost sure local convergence of G* to G°°. Indeed, on the complement 
of the event A^o, the component of (0, t/) in G°° is finite. For large enough t, it coincides 
with its component in the finite picture. Increasing t further, if necessary, we get that this 
component does not contain the oldest vertex. 
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The second statement is significantly harder to prove. We first show that on the vertex 
(0, U) is still in the infinite component of the graph \G°°, introduced in Subsection 5.2, if A 
is slightly less than 1. The remaining vertices in A!\xA! will be used at the end of the proof 
for a sprinkling argument. 

Recall that \G°° has the same law as G°° but with / replaced by A/, and thus 7 by A7. 
Taking 1 > A > we thus get \9 > 0. In particular we infer that Ac < 1 and hence, by 

Proposition]^ the mapping fi ^^0 is left-continuous at the point ^ = 1. Denote by a^oo 
the event that (0, U) is connected to infinity in the graph \G°°. Then left-continuity implies 
that P*(^oo\ A^oo) —)• 0 as A t 1- Hence it suffices to show, for fixed A G (7(37:^! 1)) that 

]P*(A^oo\^i) —^ 0. 

t^OO 

The remainder of the proof being technical and geometric, we write it for the case d = 1, for 
the sake of clarity. No argument is specific to dimension one, however, and it is not hard to 
see that the proof works, mutatis mutandis, in higher dimensions. 

We fix a small parameter a G (0,1) and a large parameter m G (l,oo). Precise constraints 
on these parameters will be given later. We show that on a^oo; the vertex (0, U) is likely to 
be connected in the finite graph aG*"* to a moderately old vertex born before time t~°‘, and 
then use the sprinkling argument to ensure that this vertex is likely to be connected to the 
core, and thus to the oldest vertex of G*™. Let \Bt be the event that neither in aG*™ nor in 
aG“ a vertex located in [—t, t] is incident to an edge of length larger than t™'/2 — t. On the 
event {U > t~°’}, which holds whp, dehne aGq, the component of (0, f7) in the subgraph of 
\G°° whose vertex set is restricted to n ([—t, t] x 1]). Let \At C {U > be the 
event that aGq is connected in \G°° by a direct edge to a vertex in [—t,t] x (0,t““). 

The proof is carried out in three steps: 

• First step: P* (\B^) —> 0. 

• Second step: P* ((a^ooH A-Bt)\ A^t) —> 0. 

>-00 

• Third step: P* ((A^tO AHt)\^tm) —)• 0. 

' t^oo 

Proof of the first step. The set xX^r\([—t, t] x (0,1]) consists of a Poisson number of vertices, 
with mean 2Xt, all with birth time uniform in (0,1). The probability that a vertex with 
birth time uniform in (0,1) is incident to an edge of length larger than K has been estimated 
in m. see Theorem 4 and its proof, and is bounded by a constant times K where rjx 
is the smallest of the three constants 1,^ — 1 and d — 1. In the robust phase, rjx £ (0,1). 
Taking K = t”^/2 — t, we easily see that whp no vertex in aT* n ([—x (0,1]) is incident 
to an edge of length larger than t"^/2 — t, as soon as the constant m is chosen larger than 
l/rjx- This reasoning works in xG°° and in aG*™ as well. This proves the first step, under 
the constraint m > l/r]x- Gl 

Proof of the second step. As P*(t/ < —)• 0, we suppose U > t~°‘. We fix 6 > 0, to be 

specified later, and split our event in two parts, depending on the value of |aGo|. 

Part (i): P*({|aGo| < t’^}r\ xAooP \Bt\ xA) —0. 

t^oo 

On aAoo, a vertex of aGq has to be connected by an edge to a vertex outside [—t, t] x 1]. 
On xBt\\At, this vertex has to be located in [— 1™/2, t™'/2]\[— t, f\. Hence one of the vertices 
of aGq is incident to an edge longer than t^~^, called a long edge. This long edge either 
links two vertices of aGq, or one vertex of aGq to a vertex located in [—t'^/2,t'^/2]\[—t,f\. 
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It is easy to see that a given vertex is unlikely to be incident to a long edge. But we can 
also prove that among all the vertices of x 1], many are incident to long edges. 

Therefore in this proof we must use the fact that xCq has few vertices (no more than t^) and 
check that these vertices are not incident to long edges. In order to do that, we explore and 
reveal the component \Co and control at each step the probability of finding a long edge. 

We first explain the exploration process and how the information about the graph is pro¬ 
gressively revealed. When exploring the neighbourhood of a vertex we use the term inedge 
to denote an edge connecting the vertex to a younger neighbour. The position and birth 
times of all vertices are all revealed at once, that is, we work conditionally on \X, and the 
remaining randomness is only encoded in the variables (V(x, x'))x,x'eAA'- For x 7^ x' G xX, 
the indegree evolution processes (.^(x, t))t>s and {Z{yd^ t))t>s' are conditionally independent. 
Start the exploration with the single vertex (0, U). Reveal all its neighbours. If (0, U) is inci¬ 
dent to a long edge, then stop the exploration and declare you found a long edge. If it is not 
the case, (0, U) is declared ‘explored’, its neighbours are declared ‘to explore’. Now choose 
a vertex left to explore. Reveal all its neighbours, except that you do not reveal whether 
it is connected by an edge to an older vertex ‘to explore’. (That edge will be checked only 
when we explore the inedges of the older vertex). If you have revealed a long edge, stop. 
Otherwise, the new neighbours you have revealed are added to the set of vertices ‘to explore’, 
and the vertex is declared explored. Continue until there are no vertices left to explore, or a 
long edge is found. An important feature of this exploration process is that it will eventually 
reveal all of xCo in at most IaCqI < steps, unless it has been stopped for finding a long 
edge. Moreover, at each step, the information gathered about the indegree evolution of the 
vertices is controlled in the following way. For each vertex that is neither explored nor to 
explore, we have revealed the absence of some inedges (those that could have linked it to an 
explored vertex, precisely). For each vertex left to explore, we have revealed the presence 
of at most one inedge, and the absence of several other inedges. Now we have to bound the 
probability of finding a long edge, conditionally on this information. 

We introduce the following notation, 

3^1 = AA'n(H,t] X (t-“,i]), 

3^2 = AA-n(([-r/ 2 ,r/ 2 ]\H,t])x(t-“,i]), 

3^3 = AA'n(([-r/2,r/2]\[-t,t])x(0,t-“)). 

Though the vertices of xCo that we explore are all in Ti, the vertices of T2 and T3 have to be 
considered as potential endvertices of long edges. With high (and even extreme) probability, 
no vertex of Ti U T2 has reached indegree more than t°‘'^ , and no vertex {y,Sy) of T3 has 
reached indegree more than Sy"^ , so we work on this event. Without extra conditioning, 
bounds on the connection probabilities are easy to establish. Indeed, if x, y G Ti U T 2 and 
e > 0, we can roughly bound the probability that they are connected by 

P(x GA y) < <^(t-“d(x,y)//(t“^)) < c,t“('+^)(^-^)d(x,y)-(^-^), (8) 

where Cg is given by the Potter bounds [6l Theorem 1.5.6]. Similarly, if x G Ti and y = 
{y,Sy) G T3, then 

P(x GA y) < (^(t-“d(x,y)//(s,-^)) < (c,f“('5-^)s-'^('5-‘^)d(x,y)-(^-^)) A 1 
< (s-^('^-^)d(x,y)-(^-=) A 1) . 


( 9 ) 
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Now suppose X G 3^1 is the vertex we are currently exploring in the exploration process, and 
y is a vertex at distance > t^~^, whose connection to x we have to check. If y G 3 ^3, then y is 
the older vertex, and its indegree evolution process is only conditioned on the nonexistence 
of some edges. In that case it only decreases the probability that it is connected to x, and 
we still can use the bound ([^. Similarly, if y G 3^i U3^2) and y is the older vertex, or x is the 
older vertex but we have not revealed the presence of any inedge of x, then the conditional 
probability of x GA y is still bounded by Q. 


We give details only for the hardest case, when x has already an inedge revealed, and x is 
older than y. We write (x, yo) for the inedge revealed, and yi, ..., y^ for the other vertices 
that have been revealed not to be linked to x. We further condition on the values of V(x, x') 
for any x' different from y, yo, yi, ..., yn, writing for the sigma-algebra generated by 
these random variables. Note that they do not determine whether x is linked to x' or not. 
However, if we know in addition that x is linked to y, yo and not to yi, ..., y^, then the 
inedges of x are all determined, as well as its indegree evolution process, which we write 
2;+(x, s)s<s^<i. If we know on the contrary that x is not linked to y, nor to yi, ..., yn, and 
only to yo, then its deterministic indegree evolution process is written 2;_(x, s)s<s^<i. The 
following computation is straightforward, 

n 

P(x GG y,x GG yo,x 4^ yi,... ,x 4^ ynlS"') =p+(y)p+(yo) ]J(1 -p+(yj)), 

n (10) 

P(x 4,^ y,x GG yo,x 4^ yi,... ,x 4^ ynl?') = (1 - P-{y))p-{yo) ]J(1 - p-(yi)), 

i=l 


where we have written p+(y) for ip{syd{x,y)/f{z+{sy—))), namely the probability that x is 
linked to y knowing that its indegree evolution process has followed until then. Similarly, 
P-{y) = ip{syd{x,y)/f{z-{sy—))). There is an easy comparison between 2;+ and Z-, namely 

2;+(x, s) = Z-(x,s) if s < Sy, 

z+(x,Sy) = z-(x,sy) + l 

Z-f-(x,s) > Z-(x,s) + l if s > Sy. 

Hence, p+(y) = p-(y), and for i G {1,... ,n}, we have (1—p_(yi)) > (1—p+(yj)). Moreover, 

P-(yo) ^ ‘f(Syod(x,yo}/f(Z-(Syo-)}) ^ piSygdix, yo)//(0)) ^ 
p+(yo) p(syod(x,yo)/f(z+(syo-))} “ (p(sygd(x,yo)/f(t‘^'y)} ~ 

where the last inequality is ensured by the Potter bounds, with a strictly positive constant 
depending on s. Now, 

F(x GG yo,x 4^ yi,... ,4^ yn|l?') 

n 

> (p+{y) + (1 -PG(y))4^““^^'^^^^)pG(yo)n(i “^’+(y*)) 

i=l 

n 

>c;t-“^('+V(yo)n(l-PG(y4)). 

i=l 


Combining with (10), we get P(x gg y|x GG yo,x 4,^ yi,...,4^ yni^') < 

Integrating with respect to the law of Vx,x' gives 


But p+{y) is always bounded by 
Fix GG ylx GG yo,x 4 ^ 


yi 


..., ^ yn) < y)-(^-^). 


( 11 ) 
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Informally, the price to pay to have a bound for the conditional probability is at most the 
multiplicative factor Adding the inequalities on every y G such that 

d(x, y) > we can bound the probability that x is incident to a long edge, conditionally 
on the beginning of the exploration process, by Ei + E 2 + E^, where 
El = Cec'-4“(27^+5-£)i-(i-fe)(5-A 137^1^ 

E 2 = ^ (|y| 

y6y2 

Es = (sy - t + A 1) . 

yeys 

This bound is independent of x. Hence the probability that the exploration process reveals 
a long edge in less than steps is bounded by t^{Ei + E 2 + E^). In other words, we have 
proved that 

F,{{\xCo\<t^}n xAoon xBt\xAt \ xX) < t\Ei + E 2 + E 3 ). 

In order to conclude (i), we have to prove that the bound is likely to be small, that is, 
goes to zero in probability. As |Ti| is whp of order t, the first term t^Ei is whp of order 
^i- 5 +£+a( 275 -i-< 5 -e)-i-fe(i+ 5 -£)^ j£ ^ ^ gj^gsen small enough, this bound goes to zero. 

For the second and third one, we show their expectation goes to zero. We have 

rl ™/2 


E 


‘ + = 2A(l-o/ {y-t + t^-y^^-Ady, 


y&y2 


which is of order ^+£). Hence E[t^ill 2 ] = s+e+a{ 2 "i&+s e)+b{5 ^ which also goes 

to zero if a, 6 and e are small enough. Finally, 


E 


ye^s 


= 2A 


rt^/2 pt “ 

/ dy ds{s-A^-A{y-t + t^yAS-e) 

Jt Jo 

^ dy ds{s-A^-Ay-(S-e) 

2A r dy(y-^/^+ C s-A^-Ay-(S-s)^g\ 

Jt^-b ^ ' 


= 2A 


< 


Finishing the calculation, this bound is 0{t with 77 = min(l /7 — 1,5 — e — 1) > 0, 

and thus E[t^Fl 3 ] = goes to zero if a and b are chosen small enough. 

Part (ii): E* ((a^oo n {UCo| > t^})\ aA*) 0. 

' t^oo 


On the event aAqo 0 {IaCIoI > we work conditionally on xCq and try to connect each 
vertex of aC'q to some vertex in [—t,t] x (0,t““]. Fix some e > 0. For any given vertex in 
A*, there exists wep a vertex of aA* with birth time in ( 0 , t““] and within distance 
because the number of such vertices follows a Poisson law of parameter 2XE. For each vertex 
of aC'q, we may then choose such a vertex, and try to connect it, with success probability 
bounded below, independently of everything else, by / f (0)). By the Potter bound, for 

each e' > 0, this is bounded below by ^ where Cg/ depends only on e'. The 

number of edges between aC'q and [—t, t] x ( 0 , t““] is therefore bounded from below by a 
binomial variable of parameters [t^\ and ^ hence it is positive whp as soon as 
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b > {6 + £'){a + e). Reducing a if necessary (as well as e and e'), we can ensure that this 
inequality is satisfied, which concludes the second step. □ 


Proof of the third step. On \At, the vertex (0, U) is connected in xG°° and within [—t, t] to a 
vertex with birth time in (0, We choose arbitrarily such a vertex x = (x, s). On xBt, all 
these connections remain in the finite graph xG^^, that is, ( 0 , C/) is connected in this finite 
graph to X. It should be enough to say x is likely to be in the 2fe-core for a well-chosen k, and 
thus connected to the oldest vertex of T^m by Proposition!^ However, due to the complex 
way we used to find x, it is not that easy to ensure it is a good vertex, or to say anything 
about its degree. This is where the sprinkling argument is used, and the reason why we have 
worked with A < 1 in the entire proof. 

We condition on aG'*™ and on the choice of x = (x, s) with s < t~°‘. The law of the graph 


U {x}, 


under this conditional law is also the (unconditioned) law of G*™(i_aT*'"u{x}, V*""), because 
the set xX^'^ is a Poisson point process of intensity 1 — A independent of aR"*"*- As a 

consequence, we know that the vertex x has whp reached degree at least /gitP) at 

time 1/2 in aG*"". As aG*"" and aG*”" are both subgraphs of G*"*, taking k > 

Proposition allows to conclude that x is whp connected to the 2A;-core and in particular 
to the oldest vertex in G*"*. Hence the same holds for (0, U), and is satisfied. □ 


6. Proof of non-robustness 


6.1. Non—robustness for 7 < 1/2. We have seen in Section 4.2 that it suffices to show 
that PG°° contains no infinite component if p is chosen small enough. We will introduce a 
notion of quick path, such that if pG“ contains an infinite component, then there exists an 
infinite quick path. Quick paths will be constructed in such a way that we can estimate their 
probability using a disjoint occurrence argument. 


First moment method based on quick paths. All the graphs we consider are locally finite, 
therefore an infinite component has to be of infinite diameter. Actually, a vertex xq of 
an infinite component is always the starting vertex of at least one infinite geodesic, that 
is an infinite path (x„)„>o, x„ = (x„,Sn), in the graph with the property that the graph 
distance between two vertices x„ and x„_|_fc is always k, for all n,k > 0. This can be 
proved with a simple diagonal argument that we leave to the reader. Note that a geodesic 
is in particular a vertex and edge self-avoiding path. Starting from any infinite geodesic 
(x„)„>o in PG°°, we now construct deterministically, in two steps, an infinite self-avoiding 
path (z,i)n>o, Zn = {zn,Un), in G°° called the quick path associated with (xn)n>o- 

First, we construct a subsequence (yn)n>o as y„ = Start with <^(0) = 0, and thus 

yo = Xq. Given (p(n), define 

Ml := {k > ip(n): 3y = (y,t) e G°° such that t > and x,^(„)Oyoxfc in G“}. 

The vertex y in the definition of A//, is called a common child of the vertices x^(^n) and x^, 
note that it is chosen in G°° and not in ^G°°. For k G AQ, the graph distance between 
the vertices x^(^n) and x^ is thus (at most) 2 in G“, while it is \k — </?(n)| in pG°°. If 
Mn is non-empty, it has to be finite and we set ip{n -|- 1 ) = maxMn- Otherwise, we set 
(p{n + 1) = ipin) -I- 1 . 




ROBUSTNESS OF SCALE-FREE SPATIAL NETWORKS 


21 


By its definition (yn)n>o satisfies the following properties: 

• for all n > 0 we have G . 

• for all n > 0 and j >2, the vertices y^ and Yn+j are not connected by an edge and 
have no common child in G°°. 

• for all n > 0, the vertices y^ and yn-i-i are either connected by an edge, or have a 
common child in G°°. 

Finally, we create a third sequence by inserting, between every pair of vertices y„ and yn-i-i 
that are not connected by an edge, the oldest common child in G°°. We obtain an infinite 
sequence (z„)„>o, which is an infinite self-avoiding path of G°°, and which we call the quick 
path associated with (x„)„>o. 

We call a vertex z„ in the quick path a regular vertex if it is older than at least one of its 
neighbours z„_i and z„+i, and we call Zji a local maximum if it does not satisfying this 
property. Similarly define the local minima. Hence a vertex z^, with n > 0, belongs to the 
sequence {yk)k>o if and only if it is regular. With this terminology the path (z„)„>o has the 
following properties: 

(i) Every regular vertex is in . The starting vertex zq is also in ^X. 

(ii) A regular vertex z„ cannot be connected by an edge to any younger vertex of the 
path, except possibly z„_i and Zn+i- 

(iii) Two regular vertices Zn and Zn+j, with n >0 and j >2, can have common children 
only if J = 2 and Zjj+i is a local maximum. In that case, z^+i is their oldest 
common child. 

Properties (ii) and (iii) depend only on the graph and not on the percolation procedure, 
and define the notion of quick paths. We also define the notion of quick paths for hnite paths, 
by restricting the quantifiers accordingly. Based on the observation that if {0,U) is in the 
infinite component of PG°° then it must be the starting point of arbitrarily long quick paths 
satisfying property (i), the first moment calculation in the next subsection shows that the 
expected number of such paths of length n goes to 0, as n t oo, if p is small enough. Note 
that, given a quick path zq-h-• • •-H-Zn in G°°, it satishes condition (i) with probability at 
most because at least ^ — 1 of the vertices on a quick path of length n are regular. 

Therefore, if we show that the expected number of quick paths of length n grows at most 
exponentially, namely it is 0{G^) for some hnite constant C, we can infer from Borel-Cantelli 
that, for any p < almost surely the component of (0, U) in PG°° is hnite. 


The expected number of quick paths has at most exponential growth. The expected number 
of quick paths of length n is given by the multiple integral 
ri 


/ IE(o,u) [|{quick paths of length n starting at (0, ii)}|] du 

Jo 

r 

duo dzi... dznlPz,o,...,z„ ((zfc)o<fc<n is a quick path in G“), 


q0,l]x(Rdx(0,l])’’ 


where zq = (0, uq) and Pzo,...,z„ is the measure P conditioned on the event {zq, ..., z„ G X}. 
Under this measure, X is simply a Poisson point process of intensity one with the points 
zo,...,Zn added. The aim is now to bound the probability of {zk)o<k<n being a quick 
path, and see how we can integrate this bound. Writing Zk = {zk,Uk) for 0 < /c < n, we 
will actually hrst integrate over (zi,..., Zn) (space integration), then over {uq, ..., Un) (time 
integration). The main step will be to prove the following proposition. 
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Proposition 16. There exists a finite constant C such that for every n > 0 and distinct 
numbers uq, ... ,Un in ( 0 , 1 ], the following inequality holds, 


dzi. ■ ■ dzn^(o,uo),{zj_,ui),...,{zr„u„){{zk,Uk))o<k<n is a Quiek path in G“) 


< 


c”n 


( 12 ) 


{Uk-i A Uk)'^{uk-i V UkY ^' 


The bound given in (12) is good in many respects. First, the proof provides a constant C 
that does not depend on the choice of the prohle function (p, but only on the attachment 
rule /. Second, the term l/{uk-\ A UkY{uk-i V UkY~^ is comparable to the probability that 
two vertices in a non-spatial equivalent of our model, with birth times Uk-i and Uk, are 
connected by an edge. But most importantly, the next lemma shows that after integrating 
over time we obtain the desired bound. 


Lemma 17. If j < 1 /2, then there exists a finite constant C such that, for any n > 0, 


duo dUn ^ UkY{Uk-l V UkY 


< C 


In 


Proof. Pick —7 > a > 7 — 1. Then carrying out the integration over Un gives 

pi rl n 

duo-■■ I dUn 




/o 


/o {uk -1 A UkY{uk-i V UkY 


<c' [ duo' 
Jo 


n —1 


dUn-l 


^n-1 


{Uk-l A UkYiuk-l V UkY ’ 


for C' = (1 + a — 7 ) ^ — (a + 7 ) The result follows from this by induction. □ 

The proposition and the lemma, combined, prove the non-robustness of (G*)t>o for any 
7 < 1/2. More precisely, if C and C are the constants given in Proposition|l 6 |and Lemma 17 


then for any p < {CC) the graph ^0°° contains no infinite component almost surely and, 
with high probability, the network YGYt>o contains no giant component. 


6.2. Proof of Proposition |16| . The proof of Proposition |16| is based on two ingredients. 
First, the definition of a quick path allows the use of a BK inequality, that splits the paths into 
small parts that interact with negative correlation. Each small part comprises no more than 
four edges, and the probability of such a path can be bounded by more or less straightforward 
integration. 


Splitting procedure. We now explain how to split the sequence zq, ... ,z„ into small parts. 
The splitting procedure depends only on their birth times uo,... ,Un, which we assume to 
be pairwise distinct, but not on the spatial positions zq, ..., Zn. The rule is simple, namely 
For i = 0,... ,n, introduce a splitting at index i if either Ui is larger 
than both Ui-\ and Ui- 2 , or m is larger than both Uj+i and Ui+ 2 - 

The boundary convention is that no condition is requested for indices outside {0,... ,n}, 
so that, for example, there is always a splitting at indices 0 and n, and there is a splitting 
at index 1 if ui > uo. We write no = 0 < ni < • • • < = n for the splitting indices in 

increasing order. These split the path into k parts. The jth part consists of the sequence 
■ ■ •,with and constituting the two boundary vertices and the other 

vertices the inside of part j. Note that the boundary vertices z„^,... belong to two 
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consecutive parts. A vertex Zj is a local maximum of a part ..., ) if rij-i < i < nj 

and we have both Ui > Ui-i (if i > Uj-i) and ut > Wj+i (if i < nj). Observe that a boundary 
vertex of a part can be a local maximum of a part without being a local maximum. We say 
a vertex Zj contributes to a part, if it belongs to, but is not a local maximum of, this part. 

Using this terminology, we observe that 

• Local maxima never contribute to any part (irrespective whether they are inside of 
a part or boundary vertices of two parts). 

• Local minima are always inside a part, and contribute to it. 

• The other vertices always contribute to exactly one part, whether they are inside it 
or at its boundary. 

For 1 < j < /s, let Aj = {z„^_^CA • • • -fAZn^.} be the event that (zn^.i, • • •, Zn^) is a path in . 
Recall that a vertex x = (x, s) is a child of Zj if Zj-H-x and Ui < s. We define Sj to be the 
(random) set of all children of vertices Zj contributing to part j, different from Zj_i and Zj+i, 
which have birth times in the interval (ui,Ui-i V tti+i). Informally, the set 'Ej contains all 
the information beyond the variables V(zj, Zj+i), Uj-i < i < nj — 1, which is needed to check 
whether Aj occurs or not. 

The following lemma justifies the splitting rule. 

Lemma 18. If (zq, ..., Zn) is a quick path, then the sets {zq, ..., z„} and Si,..., are 
pairwise disjoint. 

Proof. Property (ii) of the definition of quick paths implies that vertices of the path do 
not belong to any Ej. We now use Property (hi) and the splitting procedure to see that 
Ej and Ej/ do not intersect if j^j'. Indeed, if they intersect, this would mean that a 
vertex Zj contributing to part j and a vertex Zj/ contributing to part j' have a common 
child in {ui,Ui-i V Uj-i-i) H (rtjqUj/_i V ttj'+i). By Property (hi) we must have \ f — i\ <2. If 
\i' — i\ =2, say i' = i+ 2, then their oldest common child has to be Zj+i. As they contribute 
to different parts containing Zj+i, there must be a splitting at index i -\- 1. Hence either 
Ui-i V Uj+i = Uj-i-i, or Uj+i V Uj+3 = Uj+i. In each case, their common child with birth time 
in {ui,Ui-i V Ui+i) n (ui+ 2 ,Ui-\-i V tti+a) is older than Zj+i, and we get a contradiction. If 
— i| = 1, we can assume that f = i + 1 and there is a splitting at index i. The vertex Zj 
cannot be a local maximum, as otherwise it contributes to no part. Combining these two 
facts we get that {ui, Ui-i V Ui+i) n (ui+i, Ui V Ui-1-2) is empty and hence a contradiction. □ 

BK-inequality. We now use a version of the famous van den Berg-Kesten (BK) inequality 
to show that the probability of observing a quick path is bounded above by the product of 
the probabilities of the events Aj, for j G {1,..., A;}. The version of the BK-inequality we 
use is valid for marked Poisson processes with unit intensity on a bounded domain, see [3]. 
It states that the probability for increasing events Ei,..., to occur disjointly is bounded 
above by the product of their individual probabilities, namely 

k 

P{Ei o ... o Ek) < IfUEi). 

In the present context, an event E in the space of configurations of the marked Poisson 
processes is called increasing if, given any configuration in E, the configuration with an 
arbitrary marked point added is also in E. Disjoint oceurrence of events Ei,..., Ek is 
written as Ei o ... o Ek and defined as follows. A configuration consisting of the point set S 
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with marks {ms)s£S is in o ... o if we can find Si,... ,Sk, disjoint subsets of S, so 
that, for each j, the set Sj with marks {ms)s£Sj is in Ej. We say the marked set Sj ensures 
that Ej is realized. 


Let us see how this inequality fits in our context. Recall that we work under IPzo,...,z„ and all 
the events we consider depend only on the set of children of vertices zq, ... ,z„. Therefore, 
the events are all deterministic functionals of the following ingredients: 

(1) the random variables V(zj, Zj/), for distinct indices i,i' G {0,..., n}; 

(2) the random set X' := T’\{zo,..., z„}, which is a Poisson point process of unit in¬ 
tensity on X (0,1], together with the random marks (V(x, zq), ..., V(x, z^)) in 
[0,attached to the vertices x G X'. 

First, in order to study our problem in the framework of disjoint occurrence, we have to 
remove the dependence on the random variables V(zj,Zj/). We introduce for the 

conditional probability given V(zj, Zj/) = 0, for all \i' — f | > 2. In other words, under Pzo'',,,,,z„ 
the indegree evolution process of a vertex Zj cannot grow because of vertices Zj/ with \i' — i\ > 
2, just as if the vertex z* was not seeing them. We observe that 

]Pzo,...,z„((zn) is a quick path) < ((zn) is a quick path), 

because if (z„) is a quick path, then reducing the value of V(zj, z,/) with \i' — f| > 2 neither 
affects properties (ii) and (hi), nor does it remove edges from the quick path. We also work 
conditionally on the values V(zi, Zj+i), for z G {0,..., n — 1}. 


Second, to apply the result of [3] we need to make the underlying domain bounded. To this 
end we work with the natural finite picture approximation of the graph and of our events. 
For t finite, but large enough so that Tt contains the vertices zq, ..., z„, we construct the 
graph and we denote by the event that the jth part is a path, in the graph GL 

Now the events A\,..., A^^, are increasing events of a marked Poisson point process with unit 
intensity on x (0,1]. Applying the BK-inequality gives 




■ oAi 


V. 


Z0,Zl 


,...,Vz 


-1,Z^ 


i=i 


^Zo,Zi ; • • • ; — 


As t t oo, we know that converges locally to G° 


and thus the indicator of A) to that 


of Aj, almost surely. Moreover, similarly to Ej, we define S) as the set of all children in G^ 
of vertices z* contributing to part j, different from Zj_i and Zj+i, and with birth times in 
{ui, Ui-i V Ui+i), then we also have that, almost surely, the sets S*- coincide with the sets Ej, 
for t large enough. If zq, ..., z„ is a quick path of G°°, it is clear that for t large enough, not 
only all the A) are satisfied, but also the sets S), which ensure the events A) are satisfied. 

Consequently, the events A) have to occur disjointly for t 


are disjoint, using Lemma 18 
large enough. We get 


zo, 


(zo,...,z„ is quick I Vzo.zi, • • •, Vz„_i,2 


.) Sill’S. 

i=i 




Z0,Zl? ' 


,Vz 


1 .—l^Zn 


The event Aj depends actually only on V(z„^_^, z„^-_j+i),..., V(z„^._i, z„^.). Therefore, the 
product on the right hand side is a product of independent random variables. Taking expec¬ 
tation, and using the tower property of conditional expectation, gives 

k 

IF’zo,...,z„ (zo, • • •, Z„ is a quick path) < . 

j=i 
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Figure 3. Up to symmetry there are six types of small parts after the split¬ 
ting. Illustrated, with the index of a point on the abscissa and time on the 
ordinate, these are (i) one single edge, (ii) a V shape with two edges, (iii) a V 
shape with three edges and the end vertex of the short leg between the two 
vertices of the long leg, (iv) a V shape with three edges and both vertices of 
the long leg below the end vertex of the short leg, (v) a W shape with the 
higher end vertex on the side of the deeper valley, (vi) a W shape with the 
lower end vertex on the side of the deeper valley. 


Bound for the small parts. We have bounded the probability of observing a (long) quick 
path by the product of the probabilities of observing a path, independently for each part. 


In order to prove (12), it suffices to prove the corresponding inequality for each part j, 


dzn^-i+i ■ ■ • dzn^ Uk))o<k<n is a path in G“) 

■ n 


< 




k=l 


{Uk-i Aukp{uk-iV UkY 


Instead of treating all six possible types of parts, listed in Figure 3, we only treat the most 
complex type, numbered (v). It should become clear to the reader that we are giving bounds 
in a way that would work similarly for all the other types. To lighten notation we also 
suppose the part is the first part, that is, we suppose U3<U1<U2<R0< R4 j and show 
that 


1 


1 


I 


1 


7 1 —7 7 1—7 7 1 —7 7 1 —7 * 

U2U1 ' ' 

We introduce the canonical filtration '^t, for t £ (0,1], associated to the construction of 
the graph G°° up to time t, i.e., is the smallest cr—algebra for which the restriction of 
to vertices with birth times in (0,t] is measurable. Similarly define Observe that. 
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writing (-^z(i))te(o,i] for the indegree evolution process of vertex z, the process is 

adapted to the filtration. In the following we use ci, C2,... to denote some positive constants 
depending only on the attachment rule /. 

A change of variables from (2:1,2:2, 2:4) to (yi, y'l, ys, ys) = {z2 - 2:1,20 - ^1, ^2 - 23, 24 - 23) 

and the tower property of conditional expectation yield 

[ dzidz 2 dz 3 dz 4 ^F{Ai) 

hyihy'idy3dy3lE[E[E[E[E[E[E[lAi|5'«4-]l5'«o]l5'«o-]l5'«2]|5'«2-]l5'«i]]> 


where we have here simply written E for expectation and conditional expectation under the 
probability measure Ez°o,zi,z2,z3,z4- Rewriting the indicator of Ai as product of indicators, 
l{zo-f-)-zi-fAZ2-H-Z3}l{z3-H-Z4}, the first factor is measurable with respect to ^u^-, while the 
conditional expectation of l{z3-f-)-Z4} is equal to ‘p{u4\y'^\‘^/f{Zz3{u4—))). Using a first spatial 
integration with respect to yg, we get 

/(Z, 3 (u 4 -)) 


dy 3 E[lAi| 5 'ii 4 -] = l{zo<->zieAZ 2 f^z 3 } 


U4 


The conditional expectation of the right hand side given ^uo equals 




U4 


L 1 + Zz^(uo) 


5 , 


Mo 


< Cl l{zoOZiOZ20Z3} ^^ ^ 23 (^ 0 )) 

where the inequality follows from Lemma in the appendix. We now take conditional 
expectation given Sug-i iiote that Zz^(uo) = Zz^(uo—), and integrate in space with respect 
to y'l, to obtain the bound 

f{Zz^{uo-)). 


Cl 


uo 


- 1 {ziO-Z 2 -H>Z 3 } 


(1 + Zz 3 (uo-)) /U 4 y 


U4 


KuqJ 


The conditional expectation of this bound given can be bounded by 

C 2 ('+^-(“^))f^yi{zi»Z 2 »Z 3 )(^+^"(“^))i^V. 


Uo 


Vtt 2 / 


U4 


Vtt 2 / 


by Corollary 22 in the appendix. Further, the conditional expectation of this expression 
given Su2-^ and integrated over yi and yo, is exactly equal to 


C2 


(2 + Z^^{U2-)) /Upy f{Z^^{u2-)) /(Zz3(u2-)) (2 + Z^s{u2-)) /U4y 

Up Vu 2 / U 2 U 2 U 4 Vm2^ ’ 


and bounded by 


C 3 


(1 + ^23(112 — )) /Uo\'y (1 + .^zi('U2~)) (1 + Zz;j(u 2 — )) (1 + Zzg(u 2 — )) 


Uo 


\U 2 ^ 


U 2 


U 2 


U 4 


\U 2 J 


Using Corollary again, we bound the conditional expectation given 5 mi by 


C4- 


(^)’T| 

Vui/ U2 

^n 2 \ 

^7 (1 + Zz3(ui)) ^ 


^7 (1 + Zz3(ui)) 1 

/n4\ 

\Ul / 

U2 

Vttl / 

U 4 

\Ul / 


Finally, the expectation of that expression is bounded, using Lemma 21 again, by 

1 fUo\'^ 1 /'U 2 \'y 1 fU 2 \'^ 1 fU 4 ' 

C 5 - ' ' ' ' 


Uo \Ul/ U 2 ^Ul 


U2 VU3/ U4 \U3 


1 






















ROBUSTNESS OF SCALE-FREE SPATIAL NETWORKS 


27 


Altogether, we have proved that 

/ dzidz2dz‘idz4^'¥'{Ai) < C5 


7 1—7 7 1—7 7 1—7 7 1—7 ’ 


which gives the desired result if C is chosen at least equal to 


In this calculation, space integration is used extensively to give a simple expression for the 
density of the probability, for a vertex z with indegree k at time t—, to have a child somewhere 
with birth time in dt, 



x,(z',t) I Z{z,t-) = k) 


f{k)/t. 


It is important to perform the space integration at time t, before studying the indegree 
evolution process Z[z,s) for s < t. This method is independent of the choice of a profile 
function, showing that the argument does not involve space. But an alternative approach 
would be to use the profile function more explicitly. Because ip is regularly varying with 
index < 5 , from the Potter bounds, for any 6' < 5, there exists a finite constant c such that 
ip{x) < cx~^ , for all x > 0 . Fix a choice of such a 5 ' G (l,oo) and c > 0 . Then, for given 
z = (z, u) older than z' = (z', n'), we have 


zoz 


Z{z,u—) = k) < cf{kY ^uM(z,z^ 


-S' 


and then 

IPz.z'(z-H-z') < 1 A cE[Z(z, n'—)^ ] ^M'd(z, 

With the same outline of proof, but using Lemmaand Corollary | 22 | with p 
we could show that 

^{o,uo),izi,ui),...,{zn,u„){{{zk,Uk))o<k<n is a quick path in G°°) 
n 

< c'” 1 A ((Ufc-i A Uky{uk-1 V Ukf~'^\zk - Zk-i\)~^ , 

k=l 


5 ' 01 p = 26 ', 

( 14 ) 


for some constant c', and deduce Proposition 16 by integration over all the space variables. 
A similar bound will be used in the next subsection, without further justification. 


6 . 3 . Non—robustness for 6 > dimension one. We only need to consider 7 e 

[1/2,1). In this phase, we always have 6 > 2 . We look for ways to improve the bound from 
the previous section. Any such argument has to use the spatial structure of the network 
substantially, as the corresponding nonspatial networks are robust for 7 > 1/2. 


Let us first sketch the idea informally. Suppose that d = 1. A vertex z = (z, u) has typically 
of order u~'^ children, which may be a lot. But most of these children are typically close to z, 
namely within distance u~^, and hence their local neighbourhoods are strongly correlated. 
No matter how many vertices within distance u~^ of z belong to the component of z, it will 
not help much to connect z to vertices far away. Indeed, defining the region around z as 

Cz ■= {(•^^ u'),u > u, \z' — z\ < 2u~^ — 

see Figure 4 , we can show that the typical number of vertices outside C^, that are connected 
to z, or any other vertex in Cz, is only of order log(u“^). To estimate the probability of 
a path it therefore makes sense to consider only those edges of a quick path straddling the 
boundary of a region. This idea leads us to the notion of a trace of a quick path which we use 
to improve our bounds. Informally, forgetting about the time component and just thinking 
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about Cz as a ball around z G M, it is plausible that in dimension d = 1 few edges straddle 
the boundary of because the size of the boundary of balls in M does not grow with the 
radius. In dimension d > 1 however, if we wanted to use a similar approach, we would have 
to consider the ball of radius around a vertex. The area of its boundary is of order 

and therefore the vertices within this ball would be connected to typically at least 
vertices outside, which is already too much to carry out the proof. 


Suppose now that (0, U) belongs to an infinite component of pG°°. Then it is the starting 
point of an infinite quick path (x„)„>o in G°°, as defined in the previous section, in which 
every regular vertex is in . We define the subsequence (yn)n>o given by yn = x^(n) with 
(/9(0) = 0, and 


(p{n + l) = min{fe > (pin),Xk i 

We call {yn)o<n<m a trace of the quick path (x„)„>o. Observe that if y„ = is a 

local maximum of the quick path, then y^+i = is regular and is in . 

Therefore at least half of the vertices of the trace of a quick path are in . Arguing in 
the same way as in last subsection, we have to prove that the expected number of traces 
of length n grows at most exponentially in n. This follows from the following two results, 
which are analogous to Proposition and Lemma 17 


Proposition 19. In dimension d = 1, if \ < ^ < 1 and S > then there exists a finite 
constant G" such that, for every n > 0 and to, ■ ■ ■ ,tn G (0, 1] pairwise distinct, 


dyi . . .dynIF’yo,...,y„((yfc)o<fc<n is a trace of a quick path in G°°) 


< C"^ IT 

— Ill +1-7+7 


k=l 


-l} 1{4 > tk-l} 


(15) 


where we have written yo = (0,to) y^ = {yk,tk), and where the domain of integration is 
{{yi ,... ,yn) e K”: yk ^ for all k}. 


The improvement of this bound, compared to (12), is that if tk-i < tk, the term l/tl_it]fi^ 
has been replaced by 1/tfc- Note also that the proposition is valid only in dimension one and 
for parameters satisfying 7 G [1/2,1) and <5 > 1/(1 — 7). 


Lemma 20. There exists a finite constant G'" such that, for any n > 0, we have 


dto-- 


dtr 


n(%^ 


k=l 


+--7+7 

‘'k-Xk 


tk—l} ^ ^ f'fc—1} 


ik 


< C" 


We skip the simple proof of Lemma 20 To prove Proposition [T^ we do not bound the 


probability of a sequence being the the trace of a quick path {xk)k>o directly, but instead 
construct a third sequence {zk)k>o called an almost quick path. To this end we first define 
the enlarged region Gy^__^ around vertex y^-i = {yk-i,tk-i) by 

^Vk-i •= {(yT) e IK X (0,1]: t > tk-i,d{y,yk-i) < 2t“/^ +t“^} 

U {(y,f) G M X (0,1): f < 4_i,d(?/,yfc-i) < 


Define {zk)k>o by inserting in the infinite trace {yk)k>o a vertex y/ := x^(k)-i between the 
vertices yk-i = :>^^(k-i) and yk = x,^(fc), but only if 
• q:>{k) — 1 > ip{k — 1) and 
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Figure 4. Graphical representation of the region Cz (dark shade) and en¬ 
larged region C'^ (light shade) around the vertex z = (2,0.3) G M x (0,1], 
for 7 = 1/2. 

In other words, if is even outside the enlarged region Cy^__^ and it is not already repre¬ 
sented, we insert the vertex in connecting to y^. The infinite sequence (z^) we obtain 

is again a subsequence of the quick path (x^). Again a vertex z^ is called regular for this 
sequence if it is older than z^+i or Zfc_i, and otherwise it is called a local maximum. Observe 
that local maxima of the sequence (z^) are not necessarily local maxima of the sequence (x^), 
but regular vertices of (z^) are always regular vertices of (x^.). 

It is not hard to show that the sequence (z^) satisfies Properties (ii) and (hi) of the definition 
of quick paths. Actually, it can fail to be a quick path itself, only because it may not be 
a path, as some of the pairs (zfc,Zfc+i) are not requested to be edges of the graph. Note 
though that from the sequence (z^), one can identify the vertices y^, the inserted vertices 
y/, and which pairs (zfc,Zfc+i) are required to be edges of the graph, and which are not. A 
self-avoiding sequence satisfying Properties (ii), and (iii), such that all requested edges are 
present is called an almost quick path. 


Proof of Proposition 19. Fix n > 0 and to, ■ ■ ■ ,tn distinct times in (0, 1]. Let yi,... ,yn be 
real numbers such that, defining yo = (0,to) and y^ = {yk,tk) for 1 < A: < n, the sequence 
(yfc)o<fc<n satisfies y^ ^ Cy^_^ for all 1 < A: < n. Let 


A = {A: G {1,... ,n}: yfc ^ 

If the sequence {yk)k<n is the trace of a quick path, there must he B C A, say of cardi¬ 
nality m, and, for every k £ B, a vertex y^ G such that the sequence {'Z‘k)o<k<n+m 

obtained by inserting the vertices y/, is an almost quick path. Consequently, we have 

iP’yov.yn ((y.) o<k<n is the trace of a quick path in G°°) 

dy'fci • • • dyk^F,^o,...,z„+m ((zfc)o<fc<n-hm is an almost quick path), 

BcA^ 


where we have written ki,..., km for the ordered elements of B. The number of pairs {A, B) 
with i? C A C {1,..., n} is equal to 3”, thus in order to prove (15) it suffices to show that, 
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for any possible choice of A and B = {ki,km} C A, we have 


dyi...dyn / dy} ...dyl F; 




((zfc)o<fc<n+m is an almost quick path) 


< 


91 

3 


n 

k=l 


< 4-i} 1{4 > tk-i} 


'-fc-l'-fc 


tk 


where the domain of integration, depending on A and B, is defined by the constraints 
Yk e k ^ A, Yk ^ ^'yk-i k £ A, and y'k ^ k^ B. 

We first give a bound on the probability that (zfc)o<fc<n+m is an almost quick path in the 
same way as for a quick path. We keep now the notation {zk,Uk), resp. {yk,tk) and {y^Ak)^ 
for Zfc, resp. yk and y^, for appropriate indices k. For < 6^ < 6 we replace (14) by 

,,z„+^((zfc)o<fc<n+m is an almost quick path) 

,-S' 


■ zo, 


< 




fees (16) 

X ]))[ {itk-iAtky{tk-i'dtk9~'^\yk-yk-i\)~\ 

keA\B 

for c = c{5') some finite constant. Observe that each factor corresponds to a requested edge. 


The proof of (16) requires to split the paths into small parts and then give a bound for the 


individual probability of each part. We do not provide the detail of this proof, as this is very 
similar to the previous section. Instead, we now show how to perform the integration over 
the variables yk and y). to get an improved bound. 

We first introduce the change of variables jjk = Uk — Uk-i for 1 < A: < n and y} = u'k ~ Uk-i 
for k G B, and write yk = ivkitk) and y} = {y'k-,t'k)- The domain of integration is now 
defined by the constraints yk G 9otk i)\T*(o, 4 _i) lo^ k ^ A, yk ^ ^(04 1 ) k G A, and 
y} G C(o, 4 _j) for k G B, which is a product domain with respect to the variables jjk and 


y(.. Proposition 19 will follow if, for each k, we can integrate over yk (resp. over {y},yk) if 


k G B), the single term in the product on the right hand side of (16) involving this variable, 
and ensure the result is bounded by a constant multiple of 

1{4 < tk-i] 1{4 > tk-i] 


A-'Tf'y 




This is what we now do, considering separately the different cases. 


(A) The case k ^ A. Integrating a constant over the domain {yk‘. yk G C(0 4 i)\T*(o, 4 _i)} 
we obtain a term of order 1/tk, if tfc-i < tk, and of order l/tlr'lt'l, if tk-i > tk- 

(B) The case k G A\B. Then we have to integrate {{tk-i A tk)'^{tk-i V tk)^~'^\yk\)~^' over 
the domain {yk ■ Yk ^ C*(o 4 i)}' The reader can easily check the bound in this case. 

(Cl) The case k G B and tk-i > tk- We have to integrate IVk ~ yk\~^ 

over > 4_i, |y(,| < 2t~\ - (4)“^ and \yk\ > 2t{:}_^ + tlzltl^- Write u := \yk\ - 
and bound \yk ~ y'k\~^ by u~^ , so that the integral over gives at most a factor and 

the integral is bounded by 



ih-i) ^ 




poo 

/ u~^'du < ctl~lh^^ 


for some finite constant c. We have used that 5' > 1 and (1 — 7 )( 5 ' > 1 to obtain the right 
order for the integrals in and in u. 
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(C2) The case k ^ B and tk-i < tk- 


First, we bound the integral over younger than yfc- We have to integrate again the quantity 
Ivk - yk\~^\ but the integration is now over > tk, 
and l^fcl > Writing u := \yk\ — “^tk-i v := \yk — y'k\, 'we can similarly bound 

the integral by 



du < ct\ ^ 



v} ^'du<c'ti^^, 


for c, c' some finite constants. We have used the fact (1 — 7 )<5' > 1 to obtain the order of the 
integral in and we have used ^ > 2 to bound the integral in u. 

Second, we bound the integral over y^ older than y^. The quantity we have to integrate is now 
\yk - yk\~^' ^ the integration is over 4 -i < 4 < ^k, Wkl ^ “ ( 4 )”^! 

and \yk\ > ‘^'^k\ “b With the same notation as before we have v > and u < v, and 

the integral in u, for v fixed, gives at most a factor v, so that the integral is bounded by 




( 1 - 7 )< 5 ' 




. 1 - 5 ' 


r-ife 




for some constants c, ch We have used that 1 — <5' < —1 to obtain the right order for the 
integral in v, and that (1 — 7 ) 5 ^ — 2 > —1 to obtain the right order for the integral in 


7. Appendix: Auxiliary lemmas 


For each t G [1, 00 ] hxed, the graph G* is constructed as a growing graph with vertices placed 
in Ti and with birth times in (0,1]. The indegree of a vertex x = (x,r) at time s > r is 
denoted by Z^{x,s). The process {Z^(x,s))s>r is a time-inhomogeneous pure birth process 
started in zero at time s = r. By translation in Tj, the law of this process does not depend 
on X, and we write Z^{r,s) for Z^{{0,r), s) under the measure IP(o,r)! be., conditionally on 
the vertex we consider to be in the Poisson point process. 

This appendix provides different estimates and bounds for this process. We treat simulta¬ 
neously the cases t = 1 and t > 1 , including t = 00 , and do not stop the process at time 
s = 1. The process {Z^{r, s))s>r was already studied in [27], see in particular Lemma 8 , 
where we proved that logZ*(r, s) ~ 7 logs almost surely as s t co- Moreover, it was showr0 
in Lemma 9 that the probability of having a larger indegree than [s/r)'^ decays exponentially, 
namely 

¥{Z^{r,s) > X{s/ry) < cexp(—A/ 8 ), (17) 

for some (explicit) constant c only depending on the attachment rule /. An easy modification 
of the argument gives a similar result for the increase of the process on the interval (s, s'], i.e. 


/ l + Z\r, s') 
V 1 Z*(r, s) 


>A(s'/s)^ Z\r,s) 


< cexp(—A/ 8 ). 


(18) 


Here it is important that the bound does not depend on the value taken by Z'(r, s). A 
consequence of this exponentially decaying tail is that the moments are well-controlled, see 
the next lemma and its corollary. 


o 

This statement is actually proved under a slightly stronger assumption on that attachment rule /, namely 
that f{k) = yk + 0 ( 1 ). As explained in | 27 [ Remark 6] the results carry over to our framework at the price 
of an arbitrarily small increase of 7. This is still sufficient for our applications. 
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Lemma 21. For each p G [1, oo), there exists a constant Cp depending only on p and on the 
attachment rule, so that for every r < s < s', we have 

•1 + Z\r, s')\P 


E 


1 + Z'{r, s) 


Z\r,s) 


^ Cp ( 


P7 


Proof. For a positive random variable X, we have 

OO OO 

E[XP] <Y,^ik<X <k + l){k + ly < > k){k + ly, 

k=0 k=0 

which is bounded by an explicit finite constant if X has an explicit exponentially decaying 
bound. Apply this to the random variable 

1 + Z'{r, s') / s\'y 
1 + Z'{r, s) Vs'/ 

□ 


Corollary 22. Let X and y be vertices born before time s, and s' > s. For each p G [l,oo), 


E 


l + Z'(x,s')V /1 + Z‘(y,s') 


i + Z'(x, s)y vi + ^*(y>'S) 


Z\x,s),Z\y,s) 


^ C2p 


2py 


where c^p is as introduced in Lemma^f 


Proof. Use the Cauchy-Schwarz inequality, followed by Lemma 21 


□ 


The next lemma gives a bound on the probability of observing a small degree. 


Lemma 23. There exists a function g: (0,cx)) —)• (0,cx)) growing at infinity slower than 
polynomially, such that 

F(Z\l,r) <r''/g{r)) ^ Q. (19) 

Proof. It was proved in [27] that log(Z^(l, r))/( 7 logr) converges almost surely and in prob¬ 
ability to one. In particular, there exists a function if such that for any r] > 0 and any 
r > if{r]), we have P(Z^(l,r) < r'^U-'?)) < The function if can be chosen decreasing 
with infinite limit at zero, so that its inverse if~^ is decreasing and converging to zero. 
For any r > 1, we thus have E(Z^(l,r) < U'’)) < •^“^(r). Hence we can choose 

g{r) = which is o(r^) for any / > 0. □ 

Lemma 24. There exists a function g: (0,oo) —>■ (0,oo) growing at infinity slower than 
polynomially, such that 

sup ¥{Z\s, 1) < s~'^/g{s~y) 0. (20) 

tG(l, + oo] 

. i ^ . <s 
£ log £ — 

Proof. The supremum over t is attained when t is smallest possible, i.e. tlogt = j. Using 
that Z'{r', |) is stochastically dominated by Z'{r, ^) if r < r', we have, for t > e, 

< F{zyk^ i) < (t log t)'' / g{t log t)) 

< ¥{Z\l,k) < (t log t)'^/g{t log t)), 

where the second line follows from the scaling property. In order to prove the result, us¬ 
ing it is enough to ensure that we can choose g growing slower than polynomially such 
that (ilog£)^/c/(£logt) < {t/2)'^/g{t/2), e.g. by letting g{u) = (2 logu)"^ supi <„/2 ^ 
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We stress that the probability of having a smaller degree than expected does not decay 
exponentially. Indeed, the probability that Z°^{s, 1) = 0, i.e. the indegree of the vertex born 
at time s is still null at time 1, decays only polynomially in Hence, despite the fact 
that a vertex born at time s typically has total indegree there may well be some 

untypical vertices with much fewer inedges. 
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